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We review some aspects of weak decays of hadrons containing one heavy quark. 
The main emphasis is on B physics, in particular in the framework of the Heavy 
Quark Effective Theory. 
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1 Motivations 



In the last decades we have witnessed an impressive progress in our understand- 
ing of elementary particle physics phenomena. On the one hand, more and 
more powerful experimental tools became available, on the other a deeper in- 
sight in the theoretical understanding of elementary dynamics has been gained. 
We have reached a point where the Standard Model describing strong and elec- 
troweak interactions of elementary particles has obtained unforeseen confirma- 
tions and, simultaneously, the many questions which still remain open push us 
to look beyond it. This is the frontier of particle physics, to which all of this 
book is devoted. 

Among the many unanswered questions, the mechanism of electroweak 
symmetry breaking plays a primary role - the Higgs particle which enters 
in the Standard Model in the description of such a mechanism has not yet 
been detected. In the fermion sector, it is mostly unclear as to the observed 
hierarchy between the masses of the particles, in particular it is being debated 
whether or not we can still believe that neutrinos are massless; as far as quarks 
are concerned, problems already arise in trying to define their masses, due to 
the confining strength of QCD interactions in the long distance regime. 

Another central question is how do quarks participate in weak decays: i.e. 
the actual understanding of the mixing pattern. This last point has especially 
attracted attention both from experimentalists and theorists. The Standard 
Model describes such mixing by introducing a unitary matrix, the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix. The complex nature of this matrix is 
the only source of CP violation within the Standard Model: many efforts are 
devoted to understanding whether this is sufficient to describe all CP violating 
phenomena observed till now and all those which will hopefully be observed 
in the nearby future. This is considered one of the most fertile grounds to 
reveal physics beyond the Standard Model. The requirement of CP violating 
processes as one of the .necessary conditions to generate the observed baryon- 
antibaryon asymmetry a makes this topic extremely appealing also from the 
standpoint of cosmology. 

Weak decays are concerned with all these and many other left aside ques- 
tions. The phenomenology of such decays is impressively rich, and a huge 
number of data has been collected at the numerous experimental facilities op- 
erating in the past and at present all over the world. On the other hand, most 
of these data come from hadronic processes, and hence quantum chromody- 
namics (QCD) enters into the game. QCD is for sure the less easily tractable 
theory in elementary particle physics. While the perturbation theory has pro- 
vided us with a great predictive power in quantum electrodynamics (QED), 
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the limited realm of applicability of such an approach to the strong interac- 
tions precludes to us the possibility to give equally precise predictions in QCD. 
When facing with weak decays of hadrons, we have to take into account that 
these objects are bound states of quarks and gluons, which is a highly nonper- 
turbative effect. Therefore, the interplay of both interactions should always be 
considered, and this makes our task more difficult. 

From this point of view a particularly favorable sector is the one where 
hadrons contain a single heavy quark. The large scale is set by the heavy quark 
mass rriQ which gives us the possibility to exploit asymptotic freedom. Actu- 
ally, useful considerations can be derived when considering the limit ma. — > oo. 
This is realized within the Heavy Quark Effective Theory (HQET) □ □ which 
implements systematically new and old ideas on the dynamics of heavy quark 
systems ErQ 

The decays of heavy flavored hadrons, in particular B decays, will presum- 
ably also show a higher degree of CP violation than the kaon system, where this 
phenomenon was first observed. This is why new high luminosity machines {B 
factories) have been built to study CP violation in B decays. The possibility to 
observe a contribution of non-Standard Model particles, whatever they are, in 
loop-induced rare B decays, makes this sector appealing from both theoretical 
and experimental sides. 

These are the main motivations to study weak decays of heavy quarks, 
which we shall review in the following. The first step will be the construc- 
tion of the effective weak Hamiltonian for B decays. This is useful in many 
respects: it allows to discuss non leptonic and rare B decays; it will be a 
guide in the discussion of the Heavy Quark Effective Theory, the construction 
of which, however, follows a slightly different pattern. Next, we shall briefly 
recall the formulation of the Heavy Quark Effective Theory, focusing on the 
aspects which will be relevant for the subsequent analyses. A detailed con- 
struction of the Heavy Quark Theory is given elsewhere in this booktl As an 
application of the HQET formalism, we will discuss the relations for leptonic 
constants and semileptonic decay form factors. Very often in such cases, the 
large tuq technique has_been used in conjunction with QCD sum rulesjj which 
are reviewed elsewhere u 

The 1 / toq expansion has also proven to be a very useful tool in the com- 
putation of inclusive decays of heavy flavored hadrons. These have attracted 
considerable interest due to the possibility of reducing the hadronic uncertainty 
in the sum of exclusive modes and also because of open questions concerning 
several semileptonic B decay modes and the ratios of the beauty hadron life- 
times. We shall survey these problems in the considered context. 

Finally, we shall give some conclusions, with a short comment on some 
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unaddressed topics. 

2 Effective Weak Hamiltonian for B decays 

In the Standard Model, weak decays of quarks and leptons are mediated by 
the vector bosons W ± , Z°. The fact that these particles are very heavy with 
respect to the energies usually involved in weak decays of hadrons, makes very 
convenient the use of an effective theory in which these particles are explicitly 
integrated out of the generating functional of the theory and therefore no 
longer play a dynamical role. In particular, integrating out the W boson 
leads to the familiar pointlike four-fermion interaction originally introduced 
by Fermi in order to describe the nuclear (3 decay. The path which leads from 
the Standard Model Lagrangian to the Fermi one is nothing but an example 
of the application of the operator product expansion Ej (OPE). The virtue 
of using such a description traces back to the interplay of strong and weak 
interactions. The OPE allows us to write the amplitude relative to a certain 
weak decay process as the matrix element of an effective Hamiltonian, which 
is built as a sum of local operators weighted by suitable coefficients (Wilson 
coefficients). In this context, the operators can be viewed as effective vertices, 
and the coefficients as effective coupling constants. Since this construction 
follows the procedure of integrating out the W boson, this particle does not 
appear in the operators, but its effects are contained in the coefficients. The 
expansion can roughly be viewed as an expansion both in the inverse powers of 
M\y and in the dimension of the operators, the most important contribution 
coming from dimension six four-fermion operators, since higher contributions 
are suppressed by powers of p 2 /M^, p being a typical momentum scale of the 
considered process. In the case of the Fermi theory, it is easy to see what has 
been just stated: the local four-fermion operator appearing in the Lagrangian 
is completely independent on Mw', however, the Fermi constant contains the 
dependence on it, 

G F _ g 2 

g being the SU(2)l gauge constant. 

Most importantly, the OPE succeeds in separating the long distance regime 
from the short distance one. In fact, both the coefficients and the operators 
in the OPE depend on a scale /j,. This is defined by the fact that we integrate 
out all the particles with masses above this scale, therefore the operators de- 
scribe the physics below fx, while the coefficients take into account the physics 
above it, i.e. the short distance regime. This is why they are usually referred 
to as "short distance coefficients." However, the choice of \i is arbitrary, so 
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that the dependence on it should cancel between operators and coefficients. 
In this way we have disentangled the nonperturbative ingredients of our pro- 
cess, i.e. the matrix elements of the operators in the OPE, from the purely 
perturbative ones, the Wilson coefficients. Therefore, a good choice of fi is 
obtained requiring that the strong coupling constant is low enough to make 
meaningful the perturbative calculation. In B and D decays a common choice 
is usually fj, = 0(m,b), /i — (D(m c ), respectively. The Wilson coefficients are 
determined by requiring that, at a definite scale, the amplitudes relevant to 
a process should coincide in the full and in the effective theory, a procedure 
called matching. The perturbative calculation typically shows the appearance 
of large logarithms of the kind, 



a s ((i) In 



V m 2 



so that, even when a s (/j,) is a good expansion parameter, this product is 0(1), 
and therefore the validity of the perturbative expansion is spoiled. The power- 
ful tools of the renormalization group help us to sum up these large logarithms, 
using what is called the renormalization group improved perturbation theory. 
The leading term in this case stems from the resummation of the terms 
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(the so-called leading log approximation), 
expansion, one resums the terms 



In general, at order m in this new 



a s (tiY 



a s (n)ln 



Generally speaking, it is often insufficient to stop at the leading log approxima- 
tion, and the next-to-leading order corrections should be included, which are 
actually the truly 0(a s ) corrections with respect to the leading term and "which 
show many interesting features hidden in the leading log approximation^ 

Keeping in mind the general steps necessary to build an effective Hamil- 
tonian, let us write explicitly the effective Hamiltonian describing nonleptonic 
weak decays of B mesons: 

^ £ V Ub [dOu) Q l { b + C 2 (/x) Q™ } + h.c. (1) 



H eff 



where we have omitted penguin operators and 

Qf = {Ub)v- A [(d'u) V -A + (s'c) V -a] , 

(Uu)v-A(d'b) V -A + {Uc)v-A(s'b) v ^ A 



(2) 
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with the short notation (q% q-^V-A = <?i7ai(1 — 75)<72 and U — u, c. The primed 
fields are the weak interacting quarks, related to the mass eigenstates by the 
CKM mixing 

d! \ ( V ud V us V ub \ ( d \ 
s' = Vcd V cs V cb s . (3) 

V j \ V td V ts V tb J \ b ) 

Let us observe that QY corresponds to the same operator that we would have 
found in the Fermi pointlike approximation. In such an approximation, one 
has C% = 1, C2 — 0. QCD corrections have both the effect of introducing a 
new operator, with the same flavor content, but a different color structure, and 
of modifying the Wilson coefficients. 

We will need the general features of effective theories recalled in this sec- 
tion, in order to understand the relation between HQET and QCD, as well as 
to stress the difference arising in this situation where we just want to describe 
heavy quark dynamics and therefore cannot completely integrate them out. 



3 Heavy Quark Effective Theory 

The aim of this section is to give an intuitive picture f£ heavy quark symmetries 
and the basis of the systematic HQET formulationta 

Let us recall what is specifically meant by heavy quark. Except for the 
quark masses, QCD has a single adjustable parameter, Aqcd, which could be 
viewed as the scale separating the strong coupling regime from the weak one. 
It depends on the number of "active" flavors at a given scale. It cannot be 
precisely determined, but, roughly speaking, it can be considered as the inverse 
of the radius of a hadron, where confinement effects act to give the bound state. 
For Rhad — 1 fm, it turns out that Aqcd — 200 MeV. Therefore we can classify 
heavy quarks as those with mass larger than this parameter. Since we expect 
that typical strong interactions inside the hadron involve exchange of gluons 
with a virtuality of order Aqcd, this gives the heavy quark a special role in 
the hadron, as we shall see. 

We will therefore consider c, b, t as heavy quarks, while u, d, s will 
be considered light. However, as well as the strange quark cannot be always 
considered light in the application of chiral perturbation theory techniques, 
in the same sense the considerations which are derived within HQET should 
be applied with great caution to charm. Finally, the top quark cannot be 
considered at all, since it will decay before any bound state can be formed. 
Hence, strictly speaking, HQET is a theory of b decays. 
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In order to describe qualitatively a system containing a heavy quark Q, we 
will consider the explicit case of a heavy meson, and we will refer to its light 
antiquark plus the cloud of gluons as the light degrees of freedom. We would 
like to consider this system in the large heavy quark mass limit, itlq — > oo. In 
this limit (which is is quite close to reality) the chromomagnetic interaction 
between the heavy quark spin sq and the light degrees of freedom total angular 
momentum se dies off - it is inversely proportional to the heavy quark mass 
(just like the familiar quantum mechanical spin-orbit interaction). Therefore 
such quantities are expected to decouple in the heavy quark limit. Besides, in 
the QCD Lagrangian the dependence on the quark flavor is only contained in 
the mass term, so we also expect that the heavy quark flavor becomes irrelevant 
for rriQ — > oo. This is the intuitive content of the, by now, well-known heavy 
quark spin symmetry realized within HQET. This symmetry tells us that the 
light degrees of freedom in a heavy hadron are the same independently of the 
heavy quark spin and flavor quantum numbers. 

One notes more than once the strict analogy with atomic physics, observ- 
ing that the flavor symmetry could be compared to the fact that different 
isotopes of the same element have identical chemical properties, while the spin 
symmetry translates in atomic physics in the (almost) degeneracy of hyperfine 
levels. These observations already point out a limit of the considered frame- 
work, since they allow us to exploit the decoupling of the heavy quark in order 
to work out relations among various quantities involved in heavy hadron de- 
cays; on the other hand, they do not allow us to say anything about the light 
degrees of freedom. 

Another result that can be intuitively understood is known in the literature 
as the velocity superselection rule, which is due to Georgia The momentum of 
a heavy hadron Hq can be written as p^ H = m# while for the heavy quark: 
Pq = rriQ + fc p , where k is a residual momentum of O(Aqcd)- Since, on the 
other hand Pq = rnq Vq, one finds that, for vtlq — > oo, Vq — v 11 , i.e. strong 
interactions conserve the heavy quark velocity. 

In order to give a systematic structure to these intuitive properties, we 
would like to build up an effective field theory for large txiq. However, the 
procedure outlined in the previous paragraph is not very useful if we only 
want to describe heavy hadron interactions. To use a familiar example, if we 
were to describe processes where the W boson appears as a final or initial 
state, we could no more use the effective Fermi Lagrangian where the W is no 
more a dynamical field. For this reason, HQET is not obtained integrating out 
the heavy quark, but only the "small" components of its spinor. Therefore, we 
have to start from the QCD Lagrangian relative to the heavy quark 

C Q = Q(i p -m Q )Q (4) 
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where D is the covariant derivative in the fundamental representation, and 
redefine the heavy quark spinor Q as follows: 

Q(x) = e- im ° v - x [h v (x) + H v (x)} , (5) 

where h v and H v are defined by the application of the velocity projectors, 

K = l -^Q H v = l -^Q , (6) 

and therefore they satisfy: i>h v = h v , fl H v = —H v , Substituting Eq. (0) in 
Eq. (|J) one can derive the equation of motion for H v , 

H v = . 1 iPxhv, (7) 

IV ■ U + ZlTlQ 

where D^ = D^-v Dv^ . Equation @) can be used to eliminate H v from the 
Lagrangian, leading to the following one, equivalent to Cq: 

£ = h v iv ■ Dh v + h v i P x- „ 1 n i p ±h v . (8) 

iv ■ V + 2m,Q 

The second term can now be expanded in the powers of 1/mq, yielding 

C = h v ivDh v + -^—h v (i p ± fh v + -^—h v 9s<Ta0 n Ga \ v + O (—) . (9) 
2toq 2niQ 2 \ m Q J 

The second term represents the kinetic energy of the heavy quark due to its 
residual momentum k. The last term stems from the chromomagnetic coupling 
of the heavy quark spin to the gluon field. The fact that this term appears 
only at the order 1/mg is the origin of the spin symmetry. At the leading 
order in the heavy quark expansion, we simply have 

Chqet = h v iv ■ Dh v , (10) 

from which the following Feynman rule for the heavy quark propagator can be 
derived: 

■ <"> 

v • k 2 

as well as the one for the quark-gluon vertex, 

ig s T a Vll . (12) 

T a are the SU{3) generators, T a = A a /2, where A a stand for the Gell-Mann 
matrices. Again, we recognize that spin symmetry holds from the absence of 



gamma matrices in Eq. (p^); besides, the independence of Chqet of ttiq signals 
the flavor symmetry. Most of the results that we will report in this chapter 
exploit the symmetries of Chqet- On the other hand, the predictive power 
of the approach could be improved assessing the role of symmetry-breaking 
terms, i.e. including subleading terms in the heavy quark expansion. This can 
be done starting from the result in Eq. (Q), which gives us the possibility to 
write each QCD operator as an expansion in the powers of 1/mg. In fact, 
expanding Eq. (0), one can derive the expression for the heavy quark field Q 
in QCD at the order 1/mg, 

Q(x) = e - im ^ x (l + h v , (13) 

from which we can build any desired QCD operator, for example a current, 
at the considered order. In practice we are interested in the matrix elements 
of such operators, which we would like to express by an expansion in 1/mg. 
However, the states on which the operators act still contain a dependence on 
rtiQ, since, at the next to leading order in the heavy quark expansion, h v does 
not satisfy the same equation of motion stemming from Eq. ( JlCj ) . The way out 
is to consider once and for all Eq. (|l0|) as the HQET lagrangian; then h v exactly 
satisfies the equation jf h v — h v . The subleading terms in Eq. (|^) are treated 
as perturbations. In this way, as in the usual spirit of the perturbation theory, 
the results obtained in the presence of the perturbation can be expressed in 
terms of the unperturbated quantities. 

As a concrete example of how an asymptotic relation could be modified 
by the inclusion of I/tuq corrections, let us consider the heavy hadron masses. 
On the basis of spin symmetry, one can argue that hadrons differing only in 
the orientation of the heavy quark spin should be degenerate in mass. For 
example, this should hold in the case of the pseudoscalar and vector heavy 
flavored mesons. In the charm case the experimental values give: mjj* —mi)C^ 
140 MeV, while for the beauty we get: mj. — raj ~ 47 MeV. This is not 
surprising, since we expect heavy quark symmetry to work better in the b 
case. However, we can go further by predicting the size of the corrections: if 
these are O (I/uiq), one should have 

2 2 2 2 

m D , — m D ~ m B , — m B . 

Experimentally this is quite well confirmed, 

m 2 D , ~m 2 D ~ 0.49 GeV 2 , m|» - m| ~ 0.55 GeV 2 . 

The previous considerations allow us to express the 1 / tuq corrections in terms 
of the matrix elements of the perturbations appearing in Eq. (|l) ■ Let us define 
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the expectation values of the two subleading terms in Eq. (Q) on a hadron Hq 
as follows: 

A = 7^— (H Q \h v (i p ± ) 2 h v \H Q ) , (14) 

/4 = ^- (H Q \h v 9s<7a f a \ v \H Q ) . (15) 
ZniQ A 

Notice that the states are normalized to 2toq, so that the previous quantities 
are independent of toq . Using this notation, the mass of Hq can be written 
as follows: 

m *«=^ +a+ 4^ + °(4) ■ (l6) 

Since the chromomagnetic operator is responsible for the spin symmetry break- 
ing, and hence for the mass difference of the hadrons which differ only in the 
orientation of the heavy quark spin, it is possible to relate to such a dif- 
ference. In particular, one can write /Sq = —2 [J(J + 1) — I ] A2, where A2, as 
well as A and /ij, does not depend on toq. From the measured mass split- 
ting in the beauty mesons, where we expect our considerations to hold more 
accurately, it is possible to predict 

A 2 = -(M§, - M%) ~ 0.12 GeV 2 . 

On the other hand, in the case_pf B mesons, different controversial deter- 
minations of fj,^ exist in literature^ mainly due to the difficulty of extracting 
this parameter from experimental data. Finally, A has the role of the mass 
difference between the heavy hadron and the heavy quark in the vtiq — > 00 
limit. It is of 0{K QC d)- 

Let us observe that the relation Q(x) = e~ mQV ' x h v (x), stemming from 
Eq. ( |l3| ) in the heavy quark limit, defines the heavy quark mass rriQ. This 
turns out to be a physical parameter, given by the relation tuq = Mh q ~ A, in 
the same limit. For example, the form factors describing the decay A;, — > A c eP e 
depend on A = M\ b — m&p therefore, one could use this decay to determine 
A, so that the experimental value of M\ b would provide a determination of 
TOb. In this sense, the mass appearing in Eq. ( [l6|) can be considered as the 
nonperturbative analog of the pole mass, which, despite being a well defined 
quantity in perturbation theory, suffers from renormalon ambiguities beyond 
it El We shall see later in the discussion of inclusive decays that the adoption 
of such definition leads to the absence of I/ttiq corrections to the partonic 
prediction. 
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3. 1 Weak Decays of Heavy Mesons 

An important role is played by heavy quark spin-flavor symmetries when con- 
sidering weak decays of heavy mesons. Let us consider for example the elastic 
scattering of a pseudoscalar meson P{v) induced by an external vector current 
coupled to the heavy quark. Since in the tuq — > oo limit the light degrees of 
freedom are decoupled from the heavy quark, it looks sensible to describe this 
process as if the current were acting only on the heavy quark. We can write 

(P{v')\h v ,Yh v \P{v)) = (Q,v',s' Q \h v ^h v \Q,v,s Q ) (L,s't\L,s e ) (17) 



where L represents the light degrees of freedom. The content of Eq. (17) 
amounts to a factorization, where the last term represents the overlap of the 
light degrees of freedom. If v = v' , spin symmetry assures us that nothing 
is changed for the light degrees of freedom: the overlap is simply 5 S(jS i . This 
holds also if the heavy quark in the final state is a different one, thanks to the 
flavor symmetry. 

If v v', the decoupling of the light degrees of freedom still allows us 
to write such a factorization. However, now the light degrees of freedom in 
the final state are in general changed with respect to the initial state and the 
overlap will be a new function £(v ■ v') not depending on the spin and the 
flavor of the heavy quark, and which is normalized to 1 for v ■ v' = 1, since the 
zero component of the vector current is the generator of the flavor symmetry. 
A detailed proof of such normalization will be given later. This function is 
usually referred to as the Isgur-Wise universal form factoroEj Its universal- 
ity is apparent if we consider that it is possible to describe the semileptonic 
decays B — > D or B — > D* in terms of this single function, with the result 
of a great simplification of the theoretical study of these processes, which are 
usually described by 2 and 4 form factors, respectively. The possibility to give 
asymptotic relations between these form factors is a quite appealing feature. 
As a matter of fact, one could consider such relations as a test for models 
predicting these quantities or for the outcome of nonperturbative techniques, 
such as lattice or QCD sum rules. For example, Neubert has shown that not 
one of the commonly used quark models proposed in the literature satisfied 
the asymptotic behaviour predicted by the heavy quark symmetrytLHI 

Interestingly enough, the argument that lead us to the introduction of the 
Isgur-Wise function can be repeated in similar situations. Spin symmetry tells 
us that for each value of the light degrees of freedom total angular momentum 
se there are two degenerate states with total spin: J = sq + sg, according to 
the rules of addition of angular momentum. These ideally degenerate states 
will conveniently fit into doublets. Besides, si can be written as the sum of 
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the light antiquark spin s q plus the orbital angular momentum t of the light 
degrees of freedom with respect to the heavy quark. Therefore, a given doublet 
could be identified by the value of si and the parity of the states P = (— 1) +1 . 
If I = 0, the two degenerate states form the doublet of pseudoscalar and vector 
mesons (P, P*) (with P = D, B); these have Jf t = (0", l")i/ 2 . Therefore the 
Isgur-Wise function is the form factor describing the decay of an element of 
this doublet to another element of the same doublet. 

Let us consider the case £ = 1. Now it could be either si — 1/2 or 
S( = 3/2. The two corresponding doublets have positive parity; they are: 
Jf t = (0+,l+) 1/2 , Jf t = (l+,2+) 3/ 2. We will refer to elements of these dou- 
blets generically as (Pq, P[) and (Pi, P%), respectively. The charmed 2 + , 2 state 
has been experimentally observed and denoted as the £>2(2460) meson, with 
m D , = 2458.9 ± 2.0 MeV, T D , = 23 ± 5 MeV and m D * = 2459 ± 4 MeV, 
T D * = 25±? MeV for the neutral and charged states, respectivelyEI The HQET 
state l+ /2 can be identified with D 1 (2420), with m Dl = 2422. 2± 1.8 MeV and 
= 18.9±£-j MeV0 even though a l + y 2 component can be contained in 
such physical state due to the mixing allowed for the finite value of the charm 
quark mass. Actually, another state with J p = 1 + has been identified Ea 
with mass m D *o — (2461 ±f A ±10 ± 32) MeV. Both the states 2+ 2 and l+ 2 
decay to hadrons by d— wave transitions, which explains their narrow width; 
the strong coupling constant governing_their two-body decays can be deter- 
mined using experimental informationEj On the other hand, the strong de- 
cays of the states belonging to the sf = i + doublet (D ,D[) occur through 
s-wave transitions, with larger expected widths than in the case of the dou- 
blet | + . Indeed, analyses of the coupling constant governing the two-body 
hadronic transitions in QCD sum rules predict TLL&-^ D + ir~) ~ 180 MeV 
and T(D[° D*+ir-) ~ 165 MeVB Estimates §§ of the mixing angle a 
between D\ and D\ give a ~ 16°. Also in the case oL.CL.and B mesons 
experimental evidences of these states have been reportedE3'E3 

Again we can introduce universal form factors, one for each transition be- 
tween couples of doublets. The universal function describing the transition of 
an element of the fundamental doublet (P, P*) to the doublet (Po,P[) with 
S£ = 1/2 will be referred to as Ti/2(v ■ v'), while the transitions of the fun- 
damental doublet to the (Pi,p2) one with sg = 3/2 will be described by the 
function t 3 / 2 ( v ' v ')- However, the heavy quark symmetry does not predict 
the normalization of these functions, though a great simplification is again 
obtained, since the two r functions take the place of 14 form factors, usually 
employed to describe these transitions. 

We can stress again what has been said in the beginning: HQET is very 
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powerful in providing us with relations among different quantities entering in 
the description of heavy hadron processes, and the universal form factors are an 
important example. However, these quantities cannot be computed within the 
same context, since they represent the overlap of the light degrees of freedom. 
They are essentially nonperturbative objects, and should be computed by some 
other means. We shall see that they will play a role in the subsequent analyses. 

A very useful way to derive the structure of the matrix elements in HQET 
is to combine the members of a given doublet in a single wave function. EJ For 
example, for the fundamental doublet one can write 

H a = ^-[P;^-P a ls] , (18) 

where the first term represents the vector meson, the second the pseudoscalar 
one, and a is a light flavor index. The operators P* M , P a destroy a vector 
and a pseudoscalar meson, respectively, and contain a factor y/mp. Besides, 
v^P* — 0. Explicitly one can write 

P(v) = 

P(v,e) = V^^f~ * ■ (19) 

Finally, H a — jqHIjq. Let us now consider the transition between two heavy 
mesons P(v) — > P'(v') induced by an external current coupled to the heavy 
quarks, J = h®, Th®, where T is a product of Dirac matrices appropriate 
to the process at hand. Taking into account the previous considerations on 
the factorization of the overlap of the light degrees of freedom, it is possible 
to derive the most general decomposition for the transition matrix element, 
satisfying all the symmetry requirements (Lorentz covariance, heavy quark 
symmetry and parity), 

(P'(v')\J\P(v))=Tr[E(v,v')P'(v')rP(v)] , (20) 

where S(u, v') is a matrix containing the overlap of the light degrees of freedom 
as well as the right tensorial structure to contract free indices in the meson 
wave functions in order to reproduce the correct Lorentz structure of the matrix 
element. In the case of B — > DW semileptonic transitions, one has H(i>, v 1 ) — 
—£,(v ■ v'), and the following matrix elements can be worked out, 

{D(v')\hZ,rhl\B(v))=Z(v-v')y/Tf^E(v»+v'») , (21) 
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and 



(D*(e',v')\K,r(l-l5)h b v \B(v)) 
= i£(v ■ v') ^moniB 

x [e"^e»^ + (l + vv') e '*» - e'* v v„v"*] ■ (22) 

The two positive parity doublets introduced before can be represented analo- 
gously to Eq. (|l8|). The sg = 1/2 doublet is described by the following matrix: 



S a = ^-{P^ W -P ) , (23) 
while the se = 3/2 one is given by 

r i i 1 

(24) 



The constructions could be further extended to higher spinscJ Let us now 
demonstrate in detail the normalization of the Isgur-Wise function, since it 
will play an important role in the sequel. As it can easily be checked, for 
equal velocities, the vector current h v ^f^h v in the effective theory is conserved. 
Therefore, on account of Noether's theorem, there is a conserved charge, N = 
J d 3 xh\,{x)h v {x) . At zero recoil, using invariance under translations, one has 

d 3 x(P(v)\ht(x)h v (x)\P(v)) = f d 3 x(P{v)\ht(0)h v {0)\P{v)} . (25) 



The operator h\h v is the "number" operator, which counts the number of 
heavy quarks. Therefore, we have 

d 3 x{P(v)\hl{0)h v (0)\P{v)) = V (P{v)\P(v)) = V 2M P v° , (26) 



where V is a normalization volume. On the other hand, using Eq. (21), it 
holds also 

" d 3 x(P(v)\hl(x)h v (x)\P(v)) = V C(l) 2M P v° . (27) 



Comparing Eqs. (|2Q) - (B7J), it follows that 

= i • (28) 
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Symmetry arguments have been used in physics more than once in order to 
derive interesting phenomenological quantities. For example, the value of the 
vector coupling gy = 1 in the nuclear beta decay of the neutron, which fol- 
lows from assuming exact strong isospin invariance of nuclear interactions, has 
allowed to precisely extract V u d from such a process. Moreover, in the limit 
of exact SU(3)f symmetry, the form factor describing the semileptonic decay 
K — > irive, is normalized at zero recoil, receiving corrections only at the the 
second order in the symmetry breaking parameter m s — mu_^j This result, 
together with the analysis of such 0(m s — m u ) 2 corrections!^ has given the 
possibility to obtain an accurate determination of V us . We shall see that the 
normalization of the Isgur-Wise function will play an analogous role in the 
extraction of V c b from experimental data. 



3.2 Matching 

All our previous considerations have been derived in the heavy quark limit, 
niQ — > co. For energies much below mq we could assume that HQET rep- 
resents a good approximation to reality, in the same sense that Fermi theory 
does for energies much below M\y. However, what should be supplied in both 
cases are short distance corrections, which, as already noticed, provide new op- 
erators and modify the Wilson coefficients. In this way we reach an important 
goal: we can exploit HQET in the long distance regime, where we do not know 
how to cope with QCD; on the other hand, the corrections to the asymptotic 
predictions can be computed sistematically in perturbative QCD. 

It can be easily understood why the ultraviolet behaviour of the two the- 
ories is different: the typical terms arising in short distance computations 
involve the logarithms of the heavy quark mass, exactly as it was for the W, 
and therefore these logs diverge when taking the ttiq — > oo limit. This is par- 
ticularly important when one considers those operators, such as the vector or 
the axial- vector current, which do not require renormalization in QCD. In fact, 
these operators do require renormalization in HQET. 

Let us consider one such operator, Jqcd- We have seen that it can be ex- 
pressed in terms of operators in the effective theory, which arrange themselves 
in an expansion in the inverse wiq powers, 

Jqcd = Jhqet + — j hqet + ■ ■ ■ '> ( 29 ) 

this equality should be understood in the sense of matrix elements, since, ac- 
cording to the discussion in Sec. 0, the states on which they act are different, 
those appearing in the matrix element of the operator on the left-hand side be- 
ing dependent on itiq, on the contrary of those on which the operators on the 
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right act. Each term of the expansion receives further short distance correc- 
tions. As a consequence, at each order in l/m<g, several operators contribute 
to Eq. (p9|) , multiplied by suitable Wilson coefficients, so that each in such 
equation will become jW = ^\ Cj(fJk) Jj(fi). Let us focus our attention just 
on the first term in Eq. ((2^), and drop for simplicity the superscript and the 
subscript HQET, where from now on the operators without indices will refer 
to HQET. Equation @ becomes 

jqcd = £ cm mp) + o (^j = z r/(v) J b r e + ° (^) . 

(30) 

where we took into account the fact that HQET currents require renormaliza- 
tion by introducing a matrix of renormalization constants Zij . This is a matrix 
because in general we could have different operators with the same quantum 
numbers contributing at each order, which could mix under renormalization. 
Since both Jqcd and J bare are /x-independent quantities, the scale dependence 
should cancel between Cj(/z) and Zij(fi). This gives rise to the renormalization 
group equation for the Wilson coefficients 

m(^-7 T )c( M )=0 , (31) 

where the coefficients have been collected in a single vector C, and 7 T is the 
transpose of the matrix of anomalous dimensions, obtained by the matrix Z 
of the renormalization constants, 

1 = Z~^±Z. (32) 

At a suitable high scale m ~ 0(mQ.), one can obtain the Wilson coefficients by 
comparing the diagrams including gluon radiative corrections in the full and 
in the effective theory at the desired order in a s , since in this case there are 
no more large logs and usual perturbation theory works well. In this way, one 
can write the solution of Eq. ( |3l| ) as follows: 

C(n) = U(fx,m)C(m) . (33) 

Let us consider this solution in the case without mixing, so that 7 is a num- 
ber. This is a relevant example because, as we shall see later, this is just 
the case of the vector and axial vector currents, where no mixing occurs, and 
hence the anomalous dimension is a diagonal matrix, allowing us to apply the 



1G 



considerations which follow. In this situation, one obtains 

a s (m) - a s {n) 



U(n,m) 



a s (m) 


^0 


. «s(m) _ 





1 



4tt 



S- 



where 



7i Po ~ Pi 7o 



(34) 



(35) 



and 7o, 71 and /3o, Pi are the first two coefficients of the perturbative expansion 
of the anomalous dimension 7 and of the QCD /3-function. The approxima- 
tion of considering only the first factor in Eq. (|34|) corresponds to the leading 
log approximation, in which the large logs have been resummed. Introducing 
the one loop expression of a s and expanding again this term, one recovers the 
perturbation theory result. The second term is the next-to-leading log correc- 
tion, which requires the knowledge of the second coefficient of the anomalous 
dimension. 

If C(m) has been computed at 0(a s ), 

C(m) = C + C 1 ^, 
Air 



it is possible to write the solution in Eq. 
the heavy mass dependence, 



J33|) separating the /i dependence and 



C(ji) = C{m)K(n) 



(36) 



where 



C(m) 



1 7o/2/3o 



[a s (m)] 
[a s (fJ,)\ 



C 



1 



( C a S- 
a s {m) 

47T 



Ci) 



a s (m) 



47T 



(37) 



We shall see in the following that this formalism allows a precise study of the 
heavy hadron transitions, with relevant phcnomcnological consequences. 



4 Leptonic Decay Modes and Heavy Meson Leptonic Constants 

The leptonic decay constants are non perturbative parameters which enter 
universally in the description of bound state effects in a given particle. For a 
pseudoscalar and a vector particle, respectively, they arc defined by the matrix 
elements 

(0|A M |P(p)> - ifpp^ 
{0\V^\P*(p,X)) = ifp*mp.e M (jp,\) , 
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where e is the polarization vector of the P* meson. The leptonic constant 
fp allows the complete description of the purely leptonic decay mode of the 
pseudoscalar meson, with the decay rate given by: 



r(p-^)=^iw^(— v^(i-4y ■ (38) 

87r \mp J \ mp J 

Therefore, the leptonic decay modes are the best places where to experimen- 
tally measure such hadronic parameters; indeed, in the case of pion the con- 
stant /,r is precisely determined from the process 7r~ — > \i~Vp. However, in 
the case of heavy mesons like the B meson, things are hard. Purely leptonic 
decay channels present the difficulty of the helicity suppression factor 

m B 

that makes the purely leptonic decay mode hardly accessible in the electron and 
in the muon case (branching fractions of order 10" 10 , 10~ 7 , respectively). In 
the case of the r lepton, the helicity suppression is absent, but the experimental 
r reconstruction is a difficult task. 

One could hardly overemphasize the importance of the measurement of 
fp, due to the role played by this parameter in the description of CP violation 
in the neutral B system. Besides, as stems from Eq. (B8j), f B often appears 
together with some CKM matrix element, so that the extraction of either the 
decay constant or the relevant CKM element is affected by the uncertainty in 
the other parameter. This is especially important in the case of V u b, which is 
still affected by quite a large uncertainty. 

Adopting arguments inspired by HQET, other possibilities have been pro- 
posed, which avoid helicity suppressions!] The arguments are based on the 
quark scaling behaviour of the pseudoscalar and the vector qQ decay con- 



stants. As a matter of fact, using the trace formalism as well as Eq. dl8|) , one 



can show, in a straightforward manner, that in the limit of large b quark mass 

fp = fp* = -£= . (39) 



The constant F is another example of a universal parameter in HQET: in 
the heavy quark limit, it describes both pseudoscalar and vector heavy meson 
decays, irrespectively of their flavor. The simple relation in Eq. ( |39| ) receives 
both short distance and 0(m~^) corrections; in the leading log approximation, 
one find!S§0 

fp Imp ( a s (m c )\ 6/2 ° 



fp V m B \a s (mb) 



(40) 
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and 

fP' _ -I 2o- 8 (m q ) u-l \ 

fp ^ ' ( ' 

with P = D, B. This observation inspired the proposal @ of using radiative 
lcptonic B decays B — > ivy, in addition to purely leptonic modes, to get 
information on F and fs- As a matter of fact, B — ► processes are not 
hclicity suppressed, and their branching ratio, in the case of electron and muons 
in the final state, is estimated to be one order of magnitude larger than T(B — ► 
fiv). A determination of F by this mode, although affected by some amount 
of theoretical error, would indeed be very precious for the whole B meson 
phenomenology. 

5 Isgur-Wise Function £(y) and the Determination of \V c b\ 

Let us now consider the exclusive decays B — > D^*'£v. The weak matrix ele- 
ments relevant to these transitions can be conveniently parametrized in terms 
of form factors, 



(D{v')\V»\B{v)) 



h+(y) (v + v%+h-(y) (v-v 1 ) 



(D*(v',e)\V^\B(v)) , „ 
= 1 h v (y) e^ape v v' J 



(42) 



(D*(v> e)\A,\B(v)) 



h Al (y) (y + l)e* - [fiA 2 {y)Vf, + h As (y)v'^\ e* ■ v , 



where v and v' are the initial and final meson four-velocities, respectively, 
and y = v ■ v' . These are the form factors that, in the tuq — > 00 limit, can 
all be expressed in terms of a universal function £(?/) introduced in Sec. 2. 



Comparing Eq. fl42|) with Eqs. (|2lJ) - (g2j), one finds 

h+(y) = h v (y) = h Al (y) = h A:i (y) = £(y) , 

My) = h A2 (y) = ■ (43) 

Considering short distance corrections, all the form factors are related to a 
single function in the following, more complicated way: 

h+(y) = (c 1 (y,fi) + ^-[C 2 (y^) + C 3 (y, t i)}\ , 



M») = ^[c 2 (y,/i)-c 3 (y,/i)] z(y,»), 
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h v (y) = Cx{y,n) €(y,n) , 

h Al {y) = Cf(y,Li) £(y,fi) , 

hA 2 (y) = C$(y,n) £(y,n) , 

h A M = [cUy,ri + c!(y,»)] ■ (44) 

In the previous equations, c| 5 ^(y,/^) represent the short distance coefficients 
connecting the QCD vector (axial-vector) current to the corresponding HQET 
current; there are three such coefficients because now we have three possible 
structures, i.e. 7^, v^, v' . The /i-dependence is the same for all the functions 

, therefore, according to Eq. (|36|), one can extract such a dependence by 
writing 113 

Cf ] (y, ,i) = Cf 5 (m b , m c , y) K hh (y, /i) (45) 

where 

K hh = [a a (jM)]-**»W{l ^-Z hh (y)} , 



with dhh — §7(j/)) an( i l{y) is related to the iZElocity-dependent anomalous 
dimension of the heavy-heavy current in HQETE3 

4 

i(y) = ^[yr(y)-i], 



In 



(v + vV - i) 



Kl/) = V r^— ~ • (46) 

vr - 1 

As for the coefficient Zll. we refer the reader to the original literature where 



it has been derived. EJSOa The expressions of the coefficient functions (jf^ 
are known at the next-to-leading order£§lLHl it the leading-log approximation 
they simply read 

La s (m c )J 

with aj = — the coefficients C2 an d being zero. 

(5) 

The //-dependence of the coefficients C\ should cancel against the func- 
tion £(y,/J,), so that the form factor 

£,ren{y) = K hh (y, fi)£(y, (i) (47) 

is a renormalization group invariant function. Notice that the pattern of ra- 
diative corrections do not modify the normalization £(1) = 1. To quantify 
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this statement, let us introduce two combinations of short distance coefficients 
which will be relevant in the following discussion, 



w = c *(y = !). 1* = c !(y = !) • ( 48 ) 
i 

For the b — > c transition, the explicit calculationE! gives 

7] V = 1.025 ± 0.006, T) A = 0.986 ± 0.006 . (49) 

However, in the case of flavor conserving heavy quark currents, one has 

2as(m q ) , . 

rjv = 1, »M = 1 ^ , (50) 

which explicitly shows that the flavor conserving vector current is not renor- 
malized, assuring the £ normalization at zero recoil point. 

0(1/ mo) corrections modify this picture, introducing new operators in 
the expansion of the QCD currents in terms of the HQET ones. However, the 
matrix elements of these new operators all vanish in the zero recoil point at this 
order, so that the normalization of the matrix elements in Eq. (^2|) is further 
preserved, receiving corrections only at order 1/rriq. This is the content of 
Luke's theorem0 Such a theorem assures that the two form factors h + , /ly^in 
Eq. (^2|) receive only second order corrections in the inverse heavy quark mass 
expansion at zero recoil point, 

MD = ^ + ' 

h Al W = T]a + o(W . (51) 
\ m Qj 

However, this does not hold for all the form factors appearing in Eq. (|42|), which 
simply do not contribute at v = v' because they are multiplied by vanishing 
kinematical factorstl 

This analysis is particularly relevant since the exclusive decays B — > D^*>lv 
allow to access the CKM matrix element V c b- However, due to the explained 
pattern of corrections, only the decay to the charmed vector meson is suit- 
able, since in this case the rate, near the zero recoil point, is proportional to 
\hA 1 {y)\ 2 , which is protected by Luke's theorem. This is not the case for the 
transition to the D mcsonEa 
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The differential decay rate for the process B 
follows: 



D*iv can be written as 



^j-(B —> D*iv) = JpL 
dy 487r 3 



(m B - m D ') 2 m 3 D , \J y 1 - l(y + l) z 



(52) 



Ay m B — 2ymBmjj* + m D 



2 i 



y + 1 



(m B - m D *) 2 



\v cb \ 2 ^(y) 



where the function T 2 (y) differs from the Isgur-Wise function for short distance 
and ©(ttiq 1 ) corrections. What should be noticed is that: i) the normalization 
of T at y = 1 is known; ii) the spectrum vanishes at y = 1. The comparison 
of the experimental data with the previous result could be translated in a 
result for the product |14i,| 2 .F(l), after a suitable extrapolation to the zero 
recoil point due to the vanishing of the spectrum. This has lead to an accurate 
determination of \V c b\, which can be considered among the most important 
applications of the heavy quark symmetries. We quote the last resultpl 



\V f 



<:.b\ 



(39.8 ± 1.8 ±2.2) x 10 



-3 



(53) 



6 Determination of the Universal Form Factors from the QCD Sum 
Rules 

HQET allows one to derive model independent relations, holding in the heavy 
quark limit, mainly on the basis of the decoupling of the light degrees of free- 
dom from the spin of the heavy quark. We have noticed that in such a limit 
it is possible to introduce a whole series of nonperturbative parameters with 
a universal meaning, such as the decay constant F, or the universal form fac- 
tors describing weak decays among the members of doublets of asymptotically 
degenerate states. The usefulness of such relations is evident in the case of 
the Isgur-Wise function, where one can use symmetry arguments in order to 
predict the zero recoil point normalization and can use this information to ex- 
tract V c b from experimental data. However, we often need something more, i.e. 
the quantitative estimate of these universal quantities, which cannot be done 
within HQET, but requires the use of some nonperturbative technique, such 
as QCD sum rules or lattice QCD. In particular, the use of QCD sum rules 
in the framework of HQET revealed a very powerful tool and many HQET 
parameters have been estimated in this context. We shall review i n the fol- 
lowing subsections the determination of the Isgur-Wise function £oE3o and 
the computation of the analogous universal form factors T\/2, t 3 / 2 , describing 
the transitions of a B meson to excited positive parity charmed mesonspffi 3 ! 
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As we shall see, a preliminary part of such computations is the determina- 
tion, within the same QCD sum rule framework, of the leptonic constants of 
the heavy mesons. The functions £ and ri/2 have been computed including 
0(a s ) corrections, so that we shall have the chance to discuss the role of such 
corrections in practical situations. 

Let us summarize the basic derivation of a generic universal form factor 
using QCD sum rules. The starting point is a three-point correlator, defined 
in the effective theory, 

U lt2 (uJ,uj',y) = i 2 fdxdz e^ k ' x - kz \0\T[jf(x)J w (0)J^(zy}\0) 



where </" = qTihS, , J| = qT 2 hQ are the effective currents interpolating the 
heavy mesons with heavy quarks Q 1 and Q respectively.^] Moreover, Jy^ = 
h®, Th® is the weak current corresponding to the transition — > h®, . The 
variables k, k' are the residual momenta, obtained by the expansion of the 
heavy meson momenta in terms of the four- velocities: P = wiqv + k, P' = 
rriQiv' + k' . 

Using the analyticity of TL{u),U)' ,y) in the variables u> = 2v ■ k and u>' = 
2v' ■ k' at fixed y, one can represent the correlator ( p4] ) by a double dispersion 
relation of the form 

II( Wj *',!/)= Idvdu'- Piy ^) v) - , (55) 

J [y — u> — ie)(u' — uj — le) 

apart from possible subtraction terms. The correlator H(u>,u>',y) receives con- 
tributions from poles located at positive real values of u> and a/, corresponding 
to the physical single particle hadronic states in the spectral function p(v, z/, y). 

This contribution is proportional to the universal function S appropriate 
to the transition at hand, through the following relation: 

Tlpoie w,w ,y) cx — ■= ——= — , (56) 

(2Ai — uj — ie)(2A2 — uj' — ie) 

where /x is the renormalization scale and F\ (/x) , F2 (fi) are the effective leptonic 
constants of the heavy mesons interpolated by the currents and J| respec- 
tively, in analogy to Eq. (p9). The mass parameters Ai and A2 identify the 



a The isspei of the choice of the interpolating currents in HQET is extensively discussed by 
Dai et oD 
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position of the poles in oj and u/, and can be interpreted as binding energies 
of the heavy mesons, in accordance with Eq. Jl6|). 

The higher state contributions to p(v, v' , y) can be taken into account by 
a QCD continuum starting at some thresholds v c and v' c , and are modeled 
by the perturbative spectral function (f ert (y^v' ^y) according to the quark- 
hadron duality assumptionLj Here, pP ert is the absorptive part of the pertur- 
bative quark-triangle diagrams, with two heavy quark lines corresponding to 
the weak Q — > Q' vertex and one light quark line connecting the heavy meson 
interpolating current vertices in Eq. (|54|). At the next-to- leading order in a s , 
all possible internal gluon lines in such triangle diagrams must be considered, 
as displayed in Fig. [l| 




D 4 Da D 6 



Figure 1: Two-loop diagrams relevant for the calculation of 0(ce s ) corrections to the pertur- 
bative part of the QCD sum rule for the universal form factors. The heavy lines represent 
the heavy quark propagators in HQET. 

Therefore, for the dispersive representation in Eq. ( |55| ) in terms of hadronic 
intermediate states one assumes the ansatz 

II(u>,a/,y) =ILp i e (u!,u!',y) +U continuum (u>,u>',y) (57) 
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where, for simplicity, the dependence of the continuum contribution on the 
thresholds v c and v' c has been omitted. 

The correlator II(<jj, uj', y) can be expressed in QCD in the Euclidean re- 
gion, i.e. for large negative values of uj and uj' , in terms of perturbative and 
nonperturbative contributions 

f cf ert (v v' v) 

U(oj,u',y)= \dvdv'- 9 . ' 'JtL + ^{uj,uj',y) . (58) 

J [y — uj — ie)(y — uj — it) 

In Eq. (|8|) , Ii np represents the series of power corrections in the "small" 1 /uj 
and 1/uj' variables, determined by quark and gluon vacuum condensates or- 
dered by increasing dimension. These universal QCD parameters account for 
general properties of the nonperturbative strong interactions, for which asymp- 
totic freedom cannot be applied. The lowest dimensional ones can be obtained 
from independent theoretical sources, or fitted from other applications of QCD 
sum rules in thejpases where the hadronic dispersive contribution is particu- 
larly well-knownu In practice, since the higher dimensional condensates are 
not known, one truncates the power series and a posteriori verifies the validity 
of such an approximation. 

The QCD sum rule is finally obtained by imposing that the two repre- 
sentations of H(uj,uj' ,y), namely the QCD representation (Eq. (|58|)) and the 
pole-plus-continuum ansatz (Eq. ([57])), match in a suitable range of Euclidean 
values of uj and uj' . 

A double Borel transform in the variables uj and uj' 

-#?>=]im 7 rr(--H n > (ti->oo,lj-+ -oo,t = fixed) (59) 

t [n—ljl duj n 

(and similar for B T >) is applied to "optimize" the sum rule. As a matter of fact, 
this operation has two effects. The first one consists of factorially improving 
the convergence of the nonperturbative series, justifying the truncation proce- 
dure; the second effect enhances the role of the lowest-lying meson states while 
minimizing that of the model for the hadron continuum. The a priori undeter- 
mined mass parameters r and r' must be chosen in a suitable range of values, 
expected to be of the order of the typical hadronic mass scale (> 1 GeV), 
where the optimization is verified and, in addition, the prediction turns out 
to be reasonably stable. After the Borel transformation, possible subtraction 
terms are eliminated and Eq. (|5^) can be rewritten as 

n(r, T>,y)= f dvdv' e7 (* + ^) p(v, v' , y) . (60) 
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Equation ( p6[ ) shows that the preliminary evaluation of the constants Fi(/j,) 
and i*2(/-0 is necessary to exploit the sum rule for the determination of £. 
This is done within the same HQET QCD sum rule framework, specifically by 
evaluating two-point functions. 

In the following we shall give the results for these quantities focusing on 
the transitions B -> {D, D*) and B -> (D , D[). 



6. 1 The Isgur- Wise Function 

The main properties of the Isgur- Wise function have been already described in 
the previous sections Hprp we recall the derivation of such form factor using 
QCD sum rulesCjO'Cil As stressed above, a preliminary part of the analysis 
consists of the effective theory calculation of the leptonic constant of the heavy 
mesons belonging to the fundamental doublet (P,P*), P = D, B, since it 



enters in the sum rule for the £ function through the relation in Eq. (56). We 
have already met such a parameter in Eq. ( |39| ) ; the constant F in that equation 
is defined by 

(o\r 5 \p(y)) = , (6i) 

where J5 = qij5h v interpolates the heavy pseudoscalar meson. One considers 
therefore the two-point correlator 

= i I dW fc '*<0|T{jf (*)J|(0)}|0) . (62) 



One could have chosen equally to work with the vector current and the P* 
meson, due to the degeneracy of the two members of the doublet. 

At the next-to-leading order in renormalization group improved perturba- 
tion theory it is possible to define a renormalized scheme-independent quantity 
as follows: 

Fren = M^)] 2/9 jl - ^^[Z + 5} J F(ji) , (63) 

where, in the MS scheme, 

153-19n/ 381 - 30n/ + 287r 2 



(33 -2n/) 2 36(33- 2/1/) 

and 5 = 2/3. 

The sum rule for F(/j,) allows to evaluate also the parameter A. This can be 
done considering the logarithmic derivative i^f j;he sum rule itself with respect 
to the Borel parameter r'. The results aret 



J4fi 



0.40 ± 0.06 GeV 3/2 , A = 0.57 ± 0.07 GeV . (64) 
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The important result is that radiative corrections modify by ~ 30% the value 
of F ren , a result which is not specific of QCD sum rules, but reflects a gen- 
eral property of the considered two-point correlator, and whose origin could 
be traced back to a Coulomb interaction between quarks, since the largest en- 
hancement comes from the contribution of gluon exchange between the heavy 
and the light quark. The important role of radiativ£L.carrections in the numer- 
ical evaluation of F is confirmed by other analysescSH 

The QCD sum rule result for F can be used as an input in the three-point 
function in Eq. (|54|) in order to compute the Isgur-Wise function. In this case 
one can choose the two Dirac structures Ti, T2 to interpolate the ground state 
heavy mesons (P, P*). Therefore, the choices IY2 = ij5 (corresponding to the 
pseudoscalar meson P) or T% t 2 = 7^ — (corresponding to the vector meson 
P*) are equally acceptable, since the structure resulting from the computation 
of the trace will factor out in all the terms contributing to the sum rule. Finally, 
in Eq. (|54|), we choose J\y = K~, 7^(1 — 75)/$, in order to describe the weak 
process. 

In the lowest order in the perturbation theory, the sum rule fulfills the 
important constraint of reproducing the zero-recoil point normalization of the 
Isgur-Wise function ^ o a result which continues to hold after the inclusion 
of the a s correctioncil Hence, we observe that the QCD sum rules formulated 
in the framework of the effective theory automatically incorporate the correct 
normalization of the function £. This observation stresses the suitability of 
employing this technique in such a context. Furthermore, the inclusion of 
radiative corrections to a three-point correlator of heavy quark currents has 
been done for the first time in the HQET, where the task is simplified due 
to the simpler Feynman rules. The calculation of the loop integrals required 
to evaluate the diagrams in Fig. [l] can be done using integration by parts, 
in such a way as to reduce complicated integrals to simpler ones by recursive 
relations, and the method of differential equations developed by Kotikova The 
inclusion of radiative QCD corrections to the quark condensate had been done 
before the full inclusion of such corrections in the perturbative termE^I Also it 
has been shown cil that the calculation of the relevant loop integrals is easier 
when dealing directly with the Borel transformed sum rule. In the case of the 
Isgur-Wise function, the symmetry of the correlator allows to set equal the two 
Borel parameters: r = r', a simplification which is not allowed in the case of 
the transitions of the B meson to excited states, as we shall see in the next 
subsection. 

Let us observe that the currents employed in Eq. (jHJ) are those defined 
in the effective theory and therefore are subject to renormalization. When 
computing any universal form factor, in Eq. (p4[) there appear two heavy-light 
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currents and a heavy-heavy current; in the MS scheme the renormalization 
constants of such currents are 

Zh = 1 ~ £ ' Zhk = 1 + S 7(y) (65) 

where 

- = - +1E - ln47r, (66) 
e e 

the parameter e being defined in the D-dimensional space-time by 



D = 4 + 2e, 



reflecting the use of dimensional regularization in the employed renormaliza- 
tion scheme. j(y) has been defined in the previous section. The inclusion of 
a s corrections to the three-point correlator should reproduce the singularity 
structure of these renormalization constants. As a matter of fact, this is ex- 
plicitly verified both in the case of the £ function and for the function T\i% to 
be considered in the next section. 

We will not report here the sum rule for £, referring the reader to the 
original literaturelij What is worth noticing is that the radiative corrections 
turn out to be small, due to their cancelation in the ratio of three- to two-point 
functions. For this reason not only the normalization of £ is preserved, but 
also the impact of QCD corrections is modest for large recoil. 

In Fig. |^ we show the QCD sum rule outcome^ for the renormalized 
scheme-independent Isgur-Wise function.^] It is possible to fit the curve by 
writing 

£(j,)=£(l)[l-V (j,-l)] , (67) 

where p 2 is known as the slope of the £ function. An analp&ous decomposition 
holds for the function £ ren . Hence, the sum rule predictstHl 



p 2 ren = 0.66 ± 0.05 . (68) 

This result is in agreement with the others available in literature0El In par- 
ticular, as shown by Neubertpl it lies within the bounds which can be obtained 
on thfiJaasis of two sum rules, the Bjorken sum ruleE3 and the Voloshin sum 
ruleoE3 Such two sum rules will be discussed in the next section. 



'"The figure was taken from a review by M. Neubertl 
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Figure 2: QCD sum rule result for the Isgur-Wise function in the range of y = v ■ v' allowed 
in B — » £)(*) transitions. The dashed lines indicate the bounds on the slope at y = 1 derived 
applying the Bjorken and Voloshin sum rules. 

6.2 B decays to Excited Charmed Resonances 

It is worth analyzing other cases analogous to the determination of the Isgur- 
Wise form factor £, in particular it is interesting to calculate the universal 
form factors governing the semileptonic B meson decays into the positive 
parity charmed excited states. These higher-lying states correspond to the 
L = 1 orbital excitations in the non-relativisiic constituent quark model. Be- 
sides their theoretical relevance to HQETp3 in particular to the aspects of 
the QCD sum rule calculations, such B — ► D** semileptonic transitions (from 
now on D** will indicate the generic L — 1 charmed state) have numerous 
additional points of physical interest. Indeed, in principle these decay modes 
may account for a sizeable fraction of semileptonic i?-decays, and consequently 
they represent a well-defined set of corrections to the theoretical prediction 
that the total semileptonic B — ► X c decay rate should be saturated by the 
B — > D and B — > D* modespl in the limit mq — > oo and under the condition 
(mb — m c )/{mb + m c ) — > (the so-called small-velocity limit). Moreover, the 
shape of the inclusive differential decay distribution in the lepton energy could 
reflect contributions from the B — > D** modes. 

The issue of the contribution of such modes to the B semileptonic branch- 
ing ratio Bsl is also relevant since theoretical determinations of this quantity 
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stay above the experimental data. In fact, the parton model gives Bsl — 15%, 
with negligible 1/toq corrections.^] Perturbative QCD corrections have a more 
relevant impact, nevertheless leaving Bsl > 12.5%. Also the experimental 
results seemed controversial for a long time, since the data obtained at the 
T(AS) by the CLEO Collaboration were different fmm those obtained by the 
LEP Collaborations at the Z° peak. Recent resultsEll show an improved situ- 
ation, 

B{B^Xiv s ) = (10.63 ±0.17)% (Z° average), 
B{B^Xivi) = (10.45 ±0.21)% (T(4S) average) ; (69) 

the last number is the Particle Data Group average for the data at the T(45)0 
A possible reduction of the theoretical prediction could be obtained increasing 
the value of n c , i.e. the number of charmed hadrons produced in the B decay. 
However, since the various experimental data for these parameter show a very 
good agreement both among themselves and with theoretical predictions, it is 
likely that the solution should be found elsewhere. In particular, it is important 
to fully understand which is the contribution of the various possible decay 
modes, for example, those with D** in the final state. 

Another important result, relevant both to phenomenology and to the 
critical tests of HQET, is the relation of the B — > D** form factors at zero 
recoil to the slope p 2 of the B — ► £>W Isgur-Wise function, through the Bjorken 
sum rulep3 

n m 

In this equation n, m identify the radial excitations of states with the same J p '. 
The D** states considered here represent the lowest lying states contributing 
to the left-hand side of Eq. (frb[) . 

Of similar interest for HQET is the test of the upper bound on such uni- 
versal form factors at zero recoil,, involving the heavy meson "binding energy" 
and the D** — D mass splittingsjp'Ea 

A = 2(^^ + 2^^ A , (71) 

\ n m / 

where e k = M m - M Pq . 

Besides, the investigation of the semileptonic B transitions to excited 
charm states is an important preliminary study for the theoretical analysis 

c As we shall see in the next section, the parton model result identifies with the leading term 
in the 1/jtiq expansion for inclusive decay rates. 
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of the production of such states in nonleptonic B decaysl3't3 as well as for the 
identification of additional decay modes (such as D^^D^tt) suitable for the 
investigation of CP violating effects at B factoriesEj 

Finally, as a byproduct of the QCD sum rule calculation, theoretical pre- 
dictions for the D** meson masses can be obtained, which are obviously inter- 
esting per se. 

We have already introduced two universal form factors Ti/ 2 , t 3 / 2 describing 
the B meson transitions to the members of the two doublets with spin-parity 
J si = (0 + )l + )i/2> Jf. = (1 + , 2 + ) 3 / 2 respectively. Let us consider now the 
HQET QCD sum rule calculation of the function twj, which has been per- 
formed at the nextpto-leading order in the renormalization group improved 
perturbation theorye3 The analogous result for the .fimc.tion t 3 / 2 is available 
only at tree level in perturbative QCD corrections u3c3l£3 

The matrix elements of the semileptonic B — ► DqID and B — > D[£p tran- 
sitions can be parameterized in terms of six form factors, 



(D (v')\c^ l5 b\B(v)) 



g + (y + g-(v ~ u% . 



y/m B m Da 

(j;K, £ )|c 7M (l- 75 )6|i?W) , 
; = gv^n + e ■ v [gv 2 v^ + gv 3 v. 



fj, ' " [y v 2 -r yv 3 "fj,i , 

i gAe^ 1 e* a vl 3 v'"< , (72) 



where v and v' are four-velocities and e is the D[ polarization vector. The 
form factors gi depend on the variable y = v • v' , which is directly related to 
the momentum transfer to the lepton pair. Again, the heavy quark spin sym- 
metry allows us to relate the form factors gi(y) in Eq. ([T2|) to a single function 
Ti/2(y)E3 through short-distance coefficients, perturbatively calculable, which 
depend on the heavy quark masses mf,,m c , on y and on a mass-scale /i. The 
heavy quark spin symmetry allows us to connect the QCD vector and axial- 
vector currents to the HQET currents. At the next-to-leading logarithmic 
approximation in a s and in the infinite heavy quark mass limit, the relations 
between gi and T\i 2 are given by 

g+(y) + g-{y) = -2 (Ci 5 (j/,m) + (2/ - i)C 2 5 (y,^)) r 1/2 (y,M) , 
g+(y)-g-(y) = 2 {cf(y^) - (y - i)C|(y,^)) n /2 (y,n) , 

9Vi(y) = 2(y - 1) Ci(y,/i) Ty 2 {y,fi) , 

9v 2 {y) = -2 C 2 (y,fJ.) T 1/2 (y,fJ.) , 

gv 3 (y) = -2 (Ci(y,n) + C 3 (y,/x)J T 1/2 (y,fi) , 
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g A (y) = -2Cf(y,/i)T 1/2 (y,/i). (73) 



Following the same steps leading to Eq. (|4jj), a renormalization-group invariant 
form factor can be defined 

Tffi(y)=K hh (y,rin /2 (y,ri; (74) 

hence, in Eq. ( [73|) one can substitute C^(y,/j,) by cf\y) and T 1 / 2 {y,[i) by 
T i/2 )■ Analogous relations hold for the eight form factors parameterizing the 
matrix elements of B — ► D\tv and B — > D^lv; in this case the heavy quark 
symmetry allows to relate them to the universal function t 3 / 2 (?/)Ej The main 
difference with respect to the Isgur-Wise form factor £(y) is that one cannot 
invoke symmetry arguments to predict the normalization of both T\n{y) and 
T 3/2(y)j and therefore a calculation of the form factors in the whole kinematical 
range is required. For B — > (Do, D[)iD the physical range for the variable y is 
restricted between y — 1 and y — 1.309 — 1.326, taking into account the values 
for the mass of D , D[ (m^,^ = 2.40 - 2.45 GeV). 
Let us consider the transition B — > Dq and define 

(0\J^\P o (v))=F+(jj,) , (75) 

where JJ — qh®, represents the local interpolating current of the scalar (D ) 
meson. In analogy to Eq. ([l6|), we can write, in the heavy quark limit, 

A+ = M Do - m c . (76) 

By considering the following two-point correlator: 

*(«') = i [d^x e lfc '-(0|T[Jj'(x)Jj'(0)t]|0) , (77) 



it is possible to obtain a sum rule for the constant _F + (/i). A /i-independent 
constant F^~ en can be defined, using the relation between F + (/i) and the matrix 
element of the scalar current in full QCD, 

F+ n = Mm)] 1 W-^-Z) F+(h) (78) 

L 7T 

where Z has been defined in the previous subsection. 

Again, the sum rule for F + (^i) allows to evaluate the parameter A + , by 
considering the logarithmic derivative with respect to the Borel parameter r'. 
The following predictions are obtainedEa 

A+ = 1.0 ±0.1 GeV, F+ n = 0.7 ± 0.2 GeV^ . (79) 
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Perturbative 0(a s ) corrections represent a sizeable contributipa to the sum 
rule for F+ en , similar to the situation met in the case of fBE3c3 
Considering again the three point correlator Eq. ( |54|) and using 

J i = J v s ' = qK>, Jw = A fl = h c v ,^^ 5 h b v , and J£ = J£ = qij 5 h b v , 

it is possible to derive the sum rule for the form factor t^J£. Also in this case 
the a s correction in the perturbative term is sizeable, but it turns out to be 
partially compensated by the analogous correction in the leptonic constants 
F, F + . Notice that this is a remarkable result, not expected a priori, since 
the normalization of the form factor, for example at zero recoil, is not fixed by 
symmetry arguments. The perturbative corrections, however, do not equally 
affect the form factor for all values of the variable y, but they are sensibly y 
dependent, with the effect of increasing the slope of 7i/2jsdth respect to the 
case where they are omitted. The result is depicted in Fig. |3jpl where the curves 
refer to various choices for the continuum thresholds. The region limited by 
the curves essentially determines the theoretical accuracy of the calculation. 
Considering the y dependence, the limited range of values allowed by the mass 
difference between D and D permits the expansion near y = 1, 

r r 1 ^{y)=T 1/2 {l)(l-p\ /2 {y-l) + c l/2 {y-l) 2 ) . (80) 

A three-parameter fit to Fig. || gives 

7-1/2(1) =0.35 ±0.08 , Pi/ 2 = 2.5 ± 1.0 , c 1/2 =3±3. (81) 

An immediate application of this result concerns the prediction of the semilep- 
tonic B decay rates to -Do and D[. Using V c b — 3.9 x 10~ 2 and t{B) = 
1.56 x 10~ 12 sec, one obtains 

B(B -> D £i?) = (5±3) x 10~ 4 , B(B -> B\lv) = (7± 5) x 10~ 4 .(82) 

This means that only a very small fraction of the semileptonic B — > X c decays 
is represented by transitions into the sf — ^ + charmed doublet. 

The B meson transitions to final charmed states belonging to the si = 3/2 
doublet are described by the universal form factor T 3 / 2 (y)- This is defined by 
the following matrix elements in the heavy quark limit: 

(D 1 (v', e)|ft5,7 p (l - -f 5 )h b v \B(v)) = ^m B m Dl t 3/2 (v ■ v') 
x{[^^ e* + (£^H^ + (W-2 K] 

+i {v ""'~ 1) e liva ^{v + vT ■ (v - v'f}} , (83) 
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Figure 3: The universal form factor T^J™(y). The curves refer to choices of the threshold 

parameters: v c = 2.0 GeV, u' c = 2.5 GeV (continuous line), u c = 2.5 GeV, u' c = 3.0 GeV 
(dashed line), v c = 3.0 GeV, v' c = 3.5 GeV (dotted line). 

(D 2 (v', e)|^,7 M (l - j 5 )h b v \B(v)) = v /m B m D2 t 3/2 (v ■ v') 
x [ i ^-e^e^v^v + v'f(v - y') 7 

-V3(v ■ v> + l)e; a v a + V3e* af3 v a v%}. (84) 

The QCELsum rule analysis for this function is available in the leading or- 
der in a s c3 The procedure closely follows the one outlined above, with the 
interpolating currents 

Jr - -^«(# Y+ d v r - \fT d" 7 P )h v 

chosen to interpolate the D 2 meson with J p = 2 + , and 

JH = 975 h v 
interpolating the D\ meson with J p = 1+, 

(0\jr\D 2 (p,e)) = e^F+ 2 ^T c , 

(0|J£»|I>i(p,e)) = e^F+ 2 ^T c . (85) 
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The mass parameter analogous to those defined in Eqs. (|16[), ( |76| ) is given by 

M D2 . Dl = m c + A 3/2 . (86) 
The sum rule analysis for these quantities givesEl 
A 3/2 = 1.05 ±0.10 GeV, 



F+ 2 = 0.43 ± 0.06 GeV 5/2 . (87) 



As for the form factor r 3 / 2 (j/), the result of the sum rule is shown in Fig. [|; it 
corresponds to 

73/2(3/) = r 3 / 2 (i)[i - pt /2 {y - i)] 

with r 3 / 2 (l) ~ 0.4 and p\, 2 — 0.52. These values lead to the predictions 



B(B -> D-ylv) ~ 3.2 x 10" 



B{B -> D 2 iv) ~ 4.8 x 10" 



(89) 




Figure 4: The universal form factor T^^iy) The three curves refer to different choices of the 
threshold parameters in the perturbative contribution to the sum rule. 

Other determinations of the r functions have appeared in the literature, 
employing various versions of the constituent quark modelEil The results range 
in quite a large interval, 71/2(1) = 0.06 — 0.40, py 2 = 0.7 — 1.0 and t 3 / 2 (1) = 
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0.31 — 0.66, py 2 = 1.4 — 2.8. These results critically depend on the features of 
the models employed. 

Finally, let us mention that the universal form factors describing the B 
meson transitions to the higher doublets_with sf — | and sf = | have 
also been computed by QCD sum ruleslEl Such doublets comprise the states 
with jf = (l _ ,2~) 3 / 2 and Jf = (2~,3~) 3 / 2 , respectively. In this case, the 
application of the method to high-spin states shows several difficulties, mainly 
due to the identification of the range of the parameters needed in the sum 
rule analyses, because of the particular features of the considered states and of 
their interpolating currents. Some information could be obtained from other 
theoretical approaches, namely constituent quark models predicting the heavy 
meson spectrum. The final result, although affected by a sizeable theoretical 
uncertainty, nevertheless is useful for assessing the role of high-spin meson dou- 
blets in constituting a part of the charm inclusive semileptonic B decay width. 
Denoting the members of the considered doublets by (D*, D%), (D^, Up) cor- 
responding to sf — 3/2~ and sf — 5/2~, respectively, one finds thattj 

B(B -> Dfiut) ~ B(B D%lvi) ~ 1 x 10~ 5 . (90) 

On the other hand, the B decays to the sf = 3/2 doublet turn out to be 
negligible at the leading order in the \jm,Q expansion, due to the small value 
of the corresponding universal form factor. 

At the end of our discussion about the universal form factors describing 
B decays in the tuq — > oo limit, it should be stressed that an important issue 
concerns the role of X/ttiq corrections, which introduce subleading form factors. 
The analysis of such corrections has been carried out fqr a number of relevant 
cases, and we refer the reader to the relevant literature! 1 tttej 

7 Inclusive Heavy Hadron Decays 

We have outlined above how the use of heavy quark symmetries allows a sim- 
plified description of exclusive heavy hadron decays, both introducing universal 
form factors and allowing to classify the symmetry=hreaking contributions us- 
ing an expansion in the inverse heavy quark massLd The presence of a large 
mass also turns out to be very useful when dealing with inclusive decays, i.e. 
the transitions of a given particle into all possible final states with assigned 
quantum numbers. Such decays represent theoretically "clean" processes, since 
it is expected that in the sum over many hadrons the bound state effects, which 
are the main source of uncertainty in the calculation of exclusive transitions, 
should be eliminated by the|-averaging procedure. This expectation is based 
on the quark-hadron dualityea Besides, the possibility of computing inclusive 
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decay rates of heavy hadrons by means of an expansion in the inverse powers 
of the heavy quark mass provides in this case a reliable QCD-based systematic 
approach. 

Let us briefly summarize the main aspects of such QCD analysis of the in- 
clusive decay widths of the_hjeavy hadrons. The starting point is the transition 
operator f(Q X f -> Q)^5 

f = i J d A xT[C w {x)^ w {<S)] (91) 

describing an amplitude with the heavy quark Q having the same momentum 
in the initial and final state. Cw is the effective weak Lagrangian governing 
the decay Q — > Xf. The inclusive width of the hadron Hq can be obtained 
by averaging T over Hq and taking the imaginary part of the forward matrix 
element, 

w . I(1 = »i. (92) 

* m h q 

The reason behind calculating the right-hand side of Eq. (|9^ ) is to set up an 
operator product expansion for the transition operator T in terms of local 
operators O,, 

i 

with Oi ordered according to their dimension, and the coefficients Cj containing 
appropriate inverse powers of the heavy quark mass uiq . The lowest dimension 
operator appearing in Eq. ( |93"| ) is O3 = QQ. The next gauge and Lorentz 
invariant operator is the D — 5 chromomagnetic operator Oq — Q^a^G^Q. 

The matrix element of QQ over Hq can be obtained using the heavy quark 
equation of motion, expanded in the heavy quark mass, 

QQ = g 7 °Q + |j| - + 0(m Q 3 ) , (94) 

where 0^ is the kinetic energy operator O n = Q(iD) 2 Q. On the other hand, 
the Hq matrix element of Q~f°Q is unity (modulo the covariant normalization 
of the states) . 

The number of independent operators appearing in Eq. (|9^) increases if 
the 1/toq term is considered. Such operators are of the four-quark type 

0% = QTq qTQ (95) 
where Y is an appropriate combination of the Dirac and color matrices. 
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In this way, a complete classification of various contributions to the inclu- 
sive decay rates can be obtained for different hadrons Hq. In the expression 
for the inclusive width T(Hq — > Xf) 

T(H Q - X t ) = v{ [4 + J- + + ...] (96) 

the A{ factors which (together with I5) depend on the final state Xf, include 
perturbative short-distance coefficients and nonperturbative hadronic matrix 
elements incorporating the long range dynamics. The partonic prediction for 
the width in Eq. (|96| ) corresponds to the leading term T part (HQ —> Xf) — 
TqAq, with Aq = 1 + cJa s /n + 0(a 2 s ) and Tq oc m^; differences among the 
widths of the hadrons Hq emerge at the next to leading order in 1 /vtlq and 
are related to different values of the matrix elements of the operators Oi of 
dimension larger than three. 

It is important to notice the absence of the first order term nig 1 in Eq. (p6|), 
a result obtained by Chay, Georgi and GrinsteinE3 and Bigi, Uraltsev and 
VainshteinEll This depends on the adopted choice of the definition of the heavy 
quark mass, as explained in Sec. El. 

The occurrence of operators of the type (^5|) is an appealing feature of 
the expansion in Eq. ( p3| ) , as far as the determination of the inclusive widths 
is concerned. As a matter of fact, contrary to the D = 5 operators Oq and 
O n which are spectator blind, the D = 6 operators give different contributions 
when averaged over hadrons belonging to the same SU(3) light flavor multiplet, 
and therefore they are responsible of different lifetimes of, e.g., B~ and B s , 
Ab and 5ft. The spectator flavor dependence is related to the mechanisms of 
weak scattering and Pauli interference!^ both suppressed by the factor Wq 3 
with respect to the parton decay rate. 

In the following subsection, we shall see how the 1/mg expansion pro- 
vides a systematic framework to discuss one of the open questions in the b 
phenomenology, i.e. the problem of the ratios of beauty hadron lifetimes. 



7.1 The Problem of Beauty Hadron Lifetime Ratios 

An interesting problem of present-day heavy quark physics is represented 
by the measured difference between the Ab baryon and -EL, meson lifetimes, 
t(A°) = 1.208 ± 0.051 ps and t(B°) = 1.548 ± 0.021 psB As a matter of 
fact, the deviation from unity, at the level of 20%, of the ratio r(A{,)/r(Sd), 
r(A°)/r(-B°) = 0.780 ± 0.037^| is in contradiction with the naive expectation 

d These data were reported by W-ne Particle Data Group on the basis of the analyses reported 
by the LEP B Lifetime GroupH 
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that, at the scale of the b quark mass, the spectator model should describe 
rather accurately the decays of the hadrons containing one heavy quark. The 
ratio r(Afc)/r(i? ( i) can be computed in QCD using the previously described 
approach. We have already noticed that the first term of the expansion in 
Eq. ( |96| ) reproduces the parton model result, and therefore it contributes uni- 
versally to the inclusive decay width of all hadrons containing the same heavy 
quark. Besides, we also pointed out that the flavor of the spectator quark is 
first felt at 0(m Q 3 ). 

As for the differences in the lifetime of mesons and baryons, they could 
already arise at the level mg 2 , both due to the chromomagnetic contribution 
and to the kinetic energy term in Eq. (|9^). In particular, the kinetic energy 
term is responsible for the difference for systems where the chromomagnetic 
contribution vanishes, namely in the case of At and (the light degrees p£ 
freedom in 5— wave). However, the results of a calculation of fi^ for mesonsEj 
and baryonslZ3 support the conjecture O that the kinetic energy operator has 
practically the same matrix element when computed on such hadronic systems. 
The approximate equality of the kinetic energy operator on Bd and A;, can also 
be inferred by considering that, to the leading order in 1/mg, /i^(A(,) can be 
related Il3 to ^(B^) and to the heavy quark masses by an expression which 
assumes the charm mass m c heavy enough for a meaningful expansion in l/m c . 
Namely, 

£(A b )-ti(B d ) ~ mbmc (97) 
2(m& - m c ) 

x [(M B + 3M B « - 4M Ai ) - (M D + 3M D * - 4M Ac )] . 

Using present data and the CDF measurement M\ b = 5623 ± 5 ± 4 MevEl 
Eq. @ gives a4( a &) - (4( B d) ^ 0.002 ± 0.024 GeV 2 , where the error mainly 
comes from the error in M\ b . The QCD sum rule outcome for ^(Hq) 
is /iLLEL) ~ /i^(A(,) ~ 0.6 GeV 2 , with an estimated uncertainty of about 
30%EJ'E3 This result implies that the differences between the meson and baryon 
lifetimes should occur at the ttIq 3 level, thus involving the four-quark operators 
in Eq. (|95|). They can be classified as follows El 

Oy_ A = Qli^l qhi^Qh , 
O q s _ P = Qmh qlQr , 



2 ™ ™ " l 2 

Y qL 



t s-p = Qr — VlUl—Qr, (98) 
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1 ± 75 

with qR t L — — ^ — 1 an d the Gell-Mann matrices. 

For mesons, the vacuum saturation approximation can be used to compute 



the matrix elements of the operators in Eq. (081) , 



(B q \O v _ A \B q ) VSA = {p^)\ Bq \0%_ P \B q ) VSA = f -^^ , 

B q 

(B g \T v _ A \B g ) VSA = (B g \T$_ P \B q ) VSA = . (99) 

Therefore, the matrix elements are expressed in terms of quantities such as fs 
and the quark masses, and the resulting numerical values can be used in the 
calculation of the lifetimes, with_the only caveat concerning the accuracy of 
the factorization approximation!!.] 

The vacuum saturation approach cannot be employed for baryons; in this 
case a direct calculation of the matrix elements is required, for example using 
constituent quark models. ._. 

A simplification can be obtained Lj for Ah using color and Fierz identities 
and introducing the operators 

o%_ A = Qi^Qi <&-f<& (ioo) 

and 

&s-p = QlIr &Qr ( 101 ) 
(i and j are color indices). As a matter of fact, the A b matrix elements of 
the operators in Eq. ( p8| ) can be expressed in terms of (A b \O v _ A \A b ) and 
(Ab|0^_^|A;,), modulo 1/rriQ corrections contributing to subleading terms in 
the expression for the inclusive widths. 

The matrix element of Oy_ A and O v _ A can be parametrized as 



{U v _ A ) Ab - 2M\ b - 48 1 ' ■ 

and 



(A b \O q v _ A \A b ) = -B (A b \O q v _ A \A b ) (103) 
with B = 1 in the valence quark approximation. 



F° r /b = 200 MeV and r = 1, Eq. (102) corresponds to the value: 
(O v _ A )A b = 4.4 x 10" 3 GeV 3 . The A c matrix element of O q A has been 
computed using a bag model and a nonrelativistic quark modelO the results 
(Ov-a)a c ^ 0.75 x 10" 3 GeV 3 and (6 v _ A ) Aa ~ 2.5 x lCp 3 GeV 3 , correspond 
to r ~ CL2 and r ~ 0.6, respectively. A different model l3 is used in another 
analysislIO 
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Larger values of the matrix elements have been advocated by Rosncr 
using the values of the mass splitting — Eft and E* — E c , and assuming that 
the Aft and Eft wave functions are similar, r ~ 0.9 ±0.1, takin g — = 
M|, - M| c , or r ~ 1.8 ± 0.5 using the DELPHI measurement^ 

Information from constituent quark models have been supplemented by 
estimates based on field theoretical approaches, such as QCD sum rules and 
lattice. As a matter of fact, a large value o£r, namely r ~ 4 — 5, would explain 
the difference between r(Aft) and T(Bd)l3 The application of the QCD sum 
rule method to the calculation of the matrix element of an operator of high 
dimension presents a number of disadvantages; nevertheless, interesting and 
quite reliable information can be obtained. The resultEj is (Aft|(!?y_ JAft) ~ 
(0.4 — 1.20) x 10 -3 GeV 3 , corresponding to the parameter r in the range r ~ 
0.1 — 0.3.^ The conclusion of such analysis is that the inclusion of 1/?71q terms 
in the expression of the inclusive widths does not solve the puzzle represented 
by the difference between r(Aft) and t(Bj). As a matter of fact, using 
formulae for the lifetime ratio in the framework of the heavy quark expansion 
the QCD sum rule result together with B = 1 gives 
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T{A b )/ T (B d ) > 0.94 . (104) 

As for lattice QCD, a calculation with static b quark finds higher values for 
r, r ~ 1.2 ± 0.2, which nevertheless is not enough to explain the observed 
discrepancy in the lifetime ratioEj 

In conclusion, the expansion in the inverse heavy quark mass provides a 
systematic framework to study inclusive heavy hadron decays. In particular it 
allows to compute the corrections to the partonic picture, which identifies with 
the first term of such an expansion, in terms of the matrix elements of suitable 
operators, weighted by increasing powers of nig 1 . However, the application of 
such a formalism to the computation of the lifetime ratio r(Aft)/r(Bd) does 
not explain yet the observed ratio, leaving this issue as one of the exciting open 
problems in the phenomenology of beauty hadrons. 



8 Conclusions and Perspectives 

Weak decays of heavy quarks are especially appealing, mainly because the 
presence of a large scale, i.e. toq, allows to formulate model independent re- 
lations among nonperturbative quantities which are difficult to access in full 

e Another QCD sum rule analysisE3 finds a larger value for the parameter r. However, this 
result strongly depends on the assumption of a huge deviation from the vacuum saturation 
approximation for the four-quark condensate. 
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QCD. We have shown what we consider among the most relevant phenomeno- 
logical examples: the use of heavy quark symmetry in parametrizing weak 
matrix elements between heavy hadrons in terms of universal form factors, 
and the application of 1/jbq expansion to compute reliably inclusive decay 
rates of heavy hadrons. In both cases, the most straightforward application is 
the extraction of the CKM matrix elements V c b and V u b- The former has been 
determined with increasing accuracy, the latter, being intrinsically smaller, is 
still a subject of ongoing analyses. 

Many interesting topics have not been touched upon. The list includes 
rare B decays, among which those induced by flavor changing neutral current 
represent a fertile ground for exploring physics beyond the Standard Model; 
these are commonly viewed as those processes which can tell us something 
about new physics before the LHC era. A remarkable example are radiative 
decays induced by the b -^s-txansition, measured both in the inclusive and 
exclusive B — > K*^ modesPH which put significant constraint on the space 
of parameters of various new physics scenarios. 

Another untouched sector concerns nonlcptonic decays, where the impact 
of nonperturbative QCD effects is the highest. The effective Hamiltonian we 
have given in Sec. || is the starting point for the theoretical treatment of these 
processes. It allows us to sort out the nonperturbative quantities, the hadronic 
matrix elements of the operators in the OPE, from the perturbative ones, 
i.e. the Wilson coefficients. Having identified them, the problem of giving a 
reliable estimate of such quantities is still unsolved. Usually, either models 
or simplified approaches have been employed, among which a widely used one 
is factorization, either in-its simpler formulation, naive factorization)^ or in 
its most refined versionsEa This is quite a general problem for all nonleptonic 
decays, which is not specific of heavy hadrons. Nevertheless, the situation 
seems once more easier in this case, since it turns out that factorization works 
correctly ba^the heavy quark limit, at least for those decays usually classified 
as Class This is probably (and hopefully) a field in growth. 

Nonlcptonic decays play a central role in the analyses of CP violation in B 
decays, since the most promising channels to extract the angles of the unitarity 
triangle are just nonleptonic processes. The issue of CP violation is one of the. 
most extensively reviewed, and we refer the reader to the existing literatureEH 

Finally, I have not discussed the B s ancl_B c decays, which will play a role 
at the CERN Large Hadron Collider (LHC)E3 

It seems that in the very last years the ongoing experiments in high- 
energy physics are providing us with many exciting new results, £cpm the 
SuperKamiokande observation of the atmospheric neutrino anomalypl to the 
new measurements c3 of Re(e'/e), to the very recent direct observation of the 
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v r in July 20000 

As for B physics, the £?-factories have provided us with the very first 
results in summer 2000, thus opening the season of intense investigations in 
this sector of the elementary interactions. In a few years, the advent of LHC 
will tell us even more. The heavy hadron physics presents a number of exciting 
perspectives as well as intriguing problems still to be solved, such as the one of 
the beauty hadron lifetime ratios. Moreover, the precise measurements of CP 
violating processes and of the phases of the CKM matrix elements will provide 
us with many answers and— at the same time, will open new problems. 

In 1990, J.D. BjorkenEij pointed out that there was a great potential in the 
discovery of new symmetries in the heavy quark sector. I hope to have shown 
that he was right. 
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